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Co-Evolution of quasispecies: 
B-cell mutation rates maximize viral error catastrophes 

Christel Kamp and Stefan BornholdtB 

Institut fiir Theoretische Physik, Universitdt Kiel, Leibnizstrafle 15, D-24098 Kiel, Germany 

Co-evolution of two coupled quasispecies is studied, motivated by the competition between viral 
evolution and adapting immune response. In this co-adaptive model, besides the classical error 
catastrophe for high virus mutation rates, a second "adaptation-" catastrophe occurs, when virus 
mutation rates are too small to escape immune attack. Maximizing both regimes of viral error 
catastrophes is a possible strategy for an optimal immune response, reducing the range of allowed 
viral mutation rates to a minimum. From this requirement, one obtains constraints on B-cell muta- 
tion rates and receptor lengths, yielding an estimate of somatic hypermutation rates in the germinal 
center in accordance with observation. 
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During the past 30 years, the concept of quasispecies 
10,^ has developed into a valuable tool for modeling key 
features of molecular and viral evolution. It draws a sim- 
ple picture of sequence evolution in the presence of a 
peaked fitness function, resulting in the formation of a 
central master sequence surrounded by a cloud of mutant 
sequences, summarized as the quasispecies. A promi- 
nent feature of such systems is the occurrence of an error 
catastrophe, a sudden breakdown of stability when muta- 
tion rates get large. Recent developments of quasispecies 
models include their formulation within a statistical me- 
chanics framework and a characterization of the error 
catastrophe as a phase transition ||-|^ . 

While traditionally defined on static fitness landscapes, 
the concept of quasispecies recently has been extended to 
non-stationary fitness environments |^ . This has impor- 
tant applications, as viruses usually face quickly changing 
environments in the tight and temporal niches of their 
hosts. This new approach allows for studies of the adap- 
tive response of a quasispecies to changing external con- 
ditions. This has been studied for different choices of 
time dependent fitness functions |^,^. 

Here, we further extend this approach and study a 
virus in the environment of an adaptive immune system. 
To include the adaptive response of an immune system 
into the model, a further step has to be taken beyond a 
simple time dependent viral fitness function. An interest- 
ing observation is that motifs of immune receptors often 
form something very similar to quasispecies: A "mas- 
ter" sequence of the immune receptor, corresponding to 
the complementary epitope motif of the virus, is sur- 
rounded by a cloud of closely related receptor sequences 
that emerge from somatic hypermutation in the germi- 
nal centers |l^Jll[]. As the presence of the viral motif 
induces the proliferation of the corresponding receptor 
sequence, one can argue that this process can be repre- 
sented in the framework of quasispecies, interpreting the 
conditions that lead to proliferation of specific immune 
receptors as the fitness of the particular receptors. Let 
us, therefore, formally consider the B-cell population in- 



duced by a specific viral epitope as a quasispecies itself. 
In this paper we, therefore, study the co-evolution of two 
asymmetrically coupled quasispecies under competition. 
While the immune quasispecies is strongly attracted by 
the virus, the viral quasispecies is driven away from its 
current master sequence by the immune system. This re- 
sults in a migration through sequence space as observed 
in many infectious diseases as in HIV [|l2|,^ . In the fol- 
lowing, let us first define the model in detail. Then, dy- 
namical regimes and stability bounds are discussed which 
occur as a result of the selective forces acting on both 
sides of the system. Finally, from the perspective of an 
optimal immune response, a relationship between recep- 
tor size and mutability is derived. By maximizing the 
ranges of error catastrophes for the virus, bounds on ob- 
servable mutation rates in the immune system are de- 
rived. 

Consider a model with two quasispecies of genetic se- 
quences, one of them coding for a virus and the second 
one coding for the variable part of an immune recep- 
tor. Sequence lengths are n„ and riis, respectively, with 
bases taken from an alphabet of length A = 4. Muta- 
tion rates are quantified by the copy fidelities per base 
QviQis < 1- The time evolution of the two distributions, 
the concentration Zk{t) of viral sequence k, as well as the 
concentration of immune cells yk (t) , representing recep- 
tor coding sequence k, is described by two sets of coupled 
differential equations of the type introduced by Eigen 0] 
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Locally in sequence space, we assume a simple one-to- 
one map between viral sequence and the sequence coding 
for the immune receptor that maximally fits the viral 
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epitope (within the local neighborhood of mutated re- 
ceptors). For simplicity, both sequences share the same 
subscript / in the above formulation. Therefore, A{zi) de- 
notes the growth rate of the B-cell clone corresponding 
to receptor sequence I and depends on the concentra- 
tion of its complementary viral sequence z;. The viral 
replication rate is Bi, and its decay rate C{yi), which de- 
pends on the associated immune cell concentration. For 
viral, as well as immune receptor evolution, a transition 
probability W^i from a sequence I to sequence k by mu- 
tation is assumed, depending on the respective mutation 
rate g^, sequence length rii and the Hamming distance of 
the two sequences d(k,l). For simplicity let us assume 
nis = Uy = n with the comparable complexity of viral 
epitopes and corresponding immune receptors in mind. 
Equation (]l|) models the relative concentrations of im- 
mune receptor coding sequences, with constant overall 
population size normalized by / = -^i^i) Hi- The 
viral population (^), on the other hand, does not reach 
a constant population size, as the immune system usu- 
ally works efficiently enough for the virus not to enter 
the regime of saturation. Therefore, absolute concentra- 
tion is considered here and is the adequate quantity to 
quantify viral feed-back to immune cell proliferation. 

As viral existence in sequence space often is restricted 
to narrow niches with high fitness, let us assume a fitness 
landscape with a single peak with Bi = ^ Bm=ii = rj^ 
that moves as a result of immune system pressure. Vice 
versa, as the immune response represents a very specific 
answer to a pathogen, let us also model the immune fit- 
ness landscape to have a single peak, corresponding to 
the receptor matching the current antigenic master se- 
quence. This simplifies the co-evolution model to a few 
discrete alternatives: A{zi) = if and only if zi rep- 
resents the concentration of the viral master sequence, 
otherwise A{zi) = rjis <^ dis- Analogously, C{yi) = 5 if 
and only if yi represents the dominant immune receptor's 
concentration, otherwise C{yi) — 0. This makes the com- 
plicated couplings in the above equations tractable and 
allows us to neglect mutational backflow to the respec- 
tive master sequences. Let us first write down simplified 
equations that apply to both quasispecies (1) and (2). 
For this purpose we use non-normalized concentrations 
for quasispecies (1) also, and neglect the decay 
term in (2) for the moment. Then, each of the two qua- 
sispecies can be written in terms of the concentrations of 
a master sequence xq and of an arbitrary sequence of the 
first error class Ti xi. 

xo{t) = q^'axoit) (4) 
ii{t) = i-|g"-Vxo(t) + q'^vxiit) (5) 

The interaction between the two systems has to be spec- 
ified by extra rules on the basis of the above definition 
of growth and decay rates. To keep the model as simple 



as possible, the decay rate 5 only affects the viral se- 
quence matching the dominant immune receptor. If this 
happens to coincide with the viral master sequence, the 
viral fitness peak will effectively move. Depending on its 
strength, the former fitness peak eventually will drop be- 
low the environmental growth rate. To be specific, the 
dynamical rules of this process are defined as follows: 

1 . Once the immune system imposes a decay rate S > 
on the viral master sequence (so far stabilized at 
the viral fitness peak), the narrow niche of the virus 
is assumed to move to an arbitrary sequence of the 
first error class. 

2. The viral quasispecies adapts to the new fitness 
peak on a time scale given by the dynamical 
equations above. 

3. The fitness peak of the immune quasispecies is ad- 
justed to the new maximum of the viral distribu- 
tion. 

4. The immune system adapts to the new fitness peak 
on a second time scale determined as above. 

These steps are then iterated. While this is a strongly 
simplified picture of the co-evolutionary dynamics of two 
coupled quasispecies, the localization of the interaction 
to their respective master sequences, as well as the re- 
striction to only two growth rates for each quasispecies, 
allow a simple estimate of the dynamical regimes of the 
two coupled sets of equations of type (^,^) . Each of the 
fitness peaks is adjusted once during each cycle of dura- 
tion T = + Tis (in steps 1 resp. 3). This allows us to 
follow the arguments of Nilsson and Snoad to deter- 
mine the relative growth of the respective future master 
sequences over a full cycle t as a criterion for the quasis- 
pecies' survival 

1 XI (r) _ (e(^"^-")" - et^""-'')-) (1 - q)a 
xa{0) (A-l)(a-?7)q ' ^ ' 

This expression is applied to each one of the two qua- 
sispecies (with the respective variables) , defined over the 
full interval t — -\- Tis between two adjustments of its 
fitness peak. For a relative growth coefficient k > 1 a 
species will survive and for k < 1 it will get extinct. 

Now consider a coupled system of viral and immune 
quasispecies where the immune part exerts a selective 
pressure on the virus in the form of a non- vanishing kill or 
decay rate 5. To estimate the migration time scale of the 
virus Ti,, let us iterate the model for a full cycle r starting 
at the moment of the move of the fitness peak at f = 0. 
The relative size of the old and new master sequence 
peaks is then subsequently determined for another time 
interval t^. Let us assume a;i(0) = since the new error 
class one sequence members are mainly recruited from 
the former, weakly populated error class two. The time 
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scale Ty is given by the waiting time until the new master 
sequence population exceeds the old one: 



e(9"';.~'5)^.2;o(T) -- 



x,{t) (7) 
(eg>.-^- e^"''"")(l- qy)ay 

{X - l){av - T]y)qy 



Mutational flows between the involved sequences can be 
neglected due to the small growth of the former master 
sequence and the small size of the initial new master se- 
quence population. Assuming ay 3> rjy and gt, w 1 the 
viral adaptation timescale can be estimated to 
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qy{(Ty -■qy) + 5' 



(8) 



Similarly, for the migration time for the immune quasis- 
pecies r^g we obtain 
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(9) 



Both, Ty as well as t^s, exhibit a local minimum at spe- 
cific values of their copy fidelities q„ and g^s, mainly de- 
termined by the balance between the requirement of a 
sufficiently large initial population for the formation of 
a future master sequence and sufficiently low mutational 
losses of the new master sequence. Inserting r into the 
expressions for viral stability k„ and immune stability 
according to (^, one obtains estimates for the regimes of 
viral and immune (co-)existence. 

Fig. |l| shows these regimes in terms of the respective 
mutation rates ^ = 1 — q. The classical error catastrophe 
occurs at lower mutation rates in comparison to the static 

error threshold 1 - (^) " = 0.045 (cf. 0). This 

effect is due to additional mutational losses by migration 
and becomes large for small t. In addition, Fig. ^ shows 
different limiting behaviors for Ky and Kis for ^y and 
fiis 0, respectively, which can be summarized as 

Ky 

For the viral quasispecies one observes a second er- 
ror ("adaptability") catastrophe at small viral mutation 
rates, because a minimum viral mutation rate is needed 
to escape the decay rate S induced by the immune re- 
sponse at the viral master sequence. 
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FIG. 1. Regimes of viral and immune quasispecies 
(co-)existence, with +/— denoting stable/unstable regions of 
the respective quasispecies in dependence on mutation rates 
fiv and fiis. Parameters are cr„ — ais — 10, rjv = r^is = 1, 
S = 200, riv = riis = 50, and A = 4. A large value of S is cho- 
sen to get a good qualitative view of the system's behavior. 

The K„-surface in the /iy-/Xjs-plane, whose k = con- 
tour lines are shown in Fig. |^, is dominated by a saddle 
point: Ky{iiy) exhibits a local maximum while Ky{^iis) 
shows a local minimum. An optimal strategy for viral 
suppression is, therefore, to adjust the mutation rate fiis 
of the immune quasispecies such that k„ operates in its 
valley, with maximum regions of error catastrophes on 
both sides. One obtains the condition 
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which can be written as 



fiis -1 + Uisfiis In 



0, 



A - 1 



= 0. 



(10) 



(11) 



This mutation rate minimizes the regime of possible ex- 
istence of the viral quasispecies in Fig. |l|. Depending on 
the involved viral and immune growth rates, this range 
of allowed viral mutation rates fiy may vary (and even 
vanish for some values). 

How does this compare with experimental results? Let 
us focus on B-cells and their antibodies. Each antibody 
has at least two antigen receptors located in the variable 
regions of the antibody's heavy and light chains, each of 
which contains about 110 amino acids. Each receptor 
is coded by approximately 660 nucleotides. Antigen de- 
tection takes place in 6 subregions, the complementarity 
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determining regions (CDRs) that represent 20 — 30% of 
the antibody's variable (V-)regions ^6 17|. In the course 
of the primary immune response one observes somatic 
hypermutation in the recombined V-region genes, with 
mutational hot spots at the CDRs, resulting in an en- 
hanced affinity towards the invading antigen pO| , p^ , p^ . 
Observed mutation rates are in the range of 10~^0^^ 
mutations per base pair per generation pO|-p^. Muta- 
tion rates in the CDRs are approximately twice to tenfold 
higher than those found in the entire V-region p3| , pT[ . 
These observations are quite universal to adaptive im- 
mune systems that are common to jawed vertebrates. 
Differences are mainly found in the effectivity of selec- 
tion due to varying stages of germinal centers' expression 
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FIG. 2. Predicted mutation rates vs. receptor coding 
lengths for an optimal immune response, in comparison to 
observed rates of somatic hypermutation and observed CDR 
coding lengths. 



As Fig. H shows, the model prediction agrees well with 
the observed somatic hypermutation rates and CDR re- 
ceptor lengths. 

To summarize, the dynamics of the co-evolution of 
two coupled quasispecies has been studied. In partic- 
ular this model was formulated to provide a simple toy 
model for the co-adaptive system of viral evolution and 
immune adaptation. The model characterizes the differ- 
ent regimes of (co-)existence of viral and immune quasis- 
pecies and predicts the correct range of somatic mutation 
rates in accordance with observation. Possible extensions 
of this work are numerous, as this is only a first account 
of basic principles of co-evolving quasispecies. Analytical 
approaches beyond the simple approximation presented 
here, as well as numerical extensions may provide a more 
accurate picture of the dynamics and further possibilities 
to relate to biological data. Further applications include 
modeling HIV dynamics, e.g. by adding an overall decay 
rate representing the HIV-induced loss of CDA~^ T-cells. 



* Email address: bornholdt@theo-physik.uni-kiel.de 
[1] M. Eigen, Naturwissenschaften 58, 465 (1971). 
[2] M. Eigen and P. Schuster, The Hypercycle - A Prmct- 

ple of Natural Self-Organization (Springer Verlag, Berlin, 

1979). 

[3] M. Nowak and P. Schuster, J. Theoret. Biol. 137, 375 
(1989). 

[4] P. Tarazona, Phys. Rev. A 45, 6038 (1992). 
[5] S. Bonhoeffer and P. Stadler, J. Theoret. Biol. 164, 359 
(1993). 

[6] E. Baake, M. Baake, and H. Wagner, Phys. Rev. Lett. 

78, 559 (1997). 
[7] S. Altmeyer and J. McCaskill, Phys. Rev. Lett. 86, 5819 

(2001). 

[8] M. Nilsson and N. Snoad, Phys. Rev. Lett. 84, 191 
(2000). 

[9] C. Wilke, C. Ronnewinkel, and T. Martinetz, Phys. Rep. 
349, 395 (2001). 

[10] R. S. Harris, Q. Kong, and N. Maizels, Mutation Re- 
search 436, 157 (1999). 

[11] M. Meyer-Herrman, A. Deutsch, and M. Or-Guil, J. 
Theor. Biol. 210, 265 (2001). 

[12] S. Ganeshan, R. E. Dickover, B. Korber, Y. Bryson, and 
S. Wolinsky, J. Virol. 71, 663 (1997). 

[13] T. Allen et al.. Nature 407, 386 (2000). 

[14] B. Jones, R. Enns, and S. Rangnekar, Bull. Math. Biol. 
38, 15 (1976). 

[15] C. Thompson and J. McBride, Math. Biosci. 21, 127 
(1974). 

[16] I. Roitt, Essential Immunolgy (Blackwell Scientific Publ., 
1994). 

[17] A. Aigner and S. Neumann, Immunchemie (Fischer, 
1997). 

[18] M. Weigert, I. Cesari, S. J. Yonkovich, and M. Cohn, 

Nature 228, 1045 (1970). 
[19] G. Nossal, Cell 68, 1 (1992). 

[20] D. McKean et al., Proc. Natl. Acad. Sci. USA 81, 3180 
(1984). 

[21] C. Berek, C. Milstein, Immunol. Rev. 105, 5 (1988). 
[22] N. S. Levy, U. V. Malipiero, S. G. Lebecque, and P. J. 

Gearhart, J. Exp. Med. 169, 2007 (1989). 
[23] T. Dorner et al., J. Immunol. 158, 2779 (1997). 
[24] M. F. Diaz, M. Flajnik, Immunol. Rev. 162, 13 (1998). 
[25] M. Kasahara, Immunol. Rev. 166, 159 (1998). 
[26] M. Wilson et al., EMBO J. 11, 4337 (1992). 
[27] C.-A. Reynaud, C. Garcia, W. Hein, and J.-C. WeiU, CeU 

80, 115 (1995). 

[28] N. S. Green, M. Lin, and M. D. Scharff, Immunol. Rev. 
162, 77 (1998). 



4 



